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Abstract 

The relaxational dynamics of 1+1 dimensional directed polymer in random 

potential is studied by Monte Carlo simulations. A series of temperature 

quench experiments is performed changing waiting times tw Clear crossover 

from quasi-equilibrium behavior (t <C tw) to off-equilibrium behavior (t » t^) 

appears in the dynamical overlap function whose scaling properties are very 

similar to those found in the 3 dimensional spin-glass model. In the t ^ t^ 

part, the fluctuation dissipation theorem of the 1st kind which relates the 

response function to the tilt field with the conjugate correlation function, is 

found broken. These aging effects are brought about by the very slow growth 

of quasi-equilibrium domain driven by successive loop- excitations of various 

sizes, which form complex network structures. 
02.50.Ey, 74.60.Ge, 75.10.Nr. 
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I. INTRODUCTION 



Directed polymer in random media (DPRM) [0 is one of the simplest statistical me- 
chanical models in which quenched randomness plays non-trivial roles as in spin-glasses. It 
is an effective model of an elastic string in random environment, such as a vortex line in 
oxide cuprate superconductor which penetrates the stacked Cu02 layers with point defects 
scattered randomly over the layers. In transverse dimensions d less than two, there is no free 
phase but a pinned phase, i.e. the polymer is mostly pinned around the ground state and 
cannot move freely at any finite temperature. However there are anomalously large ther- 
mal fluctuation due to the thermal hoppings between the excited states which are nearly 
degenerate with the ground state but located far away. They are analogous to the droplet 
excitations in spin-glass phase and bring about non-trivial effects in the pinned phase 
0-^. One naturally expects that they have dramatic effects also on the dynamics. Actually, 
it has been argued in the theory of the transport problem in the vortex glass phase that 
they are responsible for the nonlinear current- voltage response 0]. 

Our interest in the present study is the slow relaxation to equilibrium due to such anoma- 
lous thermal hoppings. We expect that they bring about aging effects, which are to some 
extent similar to those found in the spin-glass phase. In order to clarify this possibility, we 
performed temperature quench experiments by a simple heatbath Monte Carlo dynamics. 
The procedure mimic the so called IRM experiments in spin-glasses [0 but the system is 
perturbed by a small tilt field which drives one end of the polymer instead of magnetic field 
in spin-glasses which drives the whole spins . 

We have found clear evidences of the aging effects. One is the systematic waiting time de- 
pendence on the dynamical overlap function. The crossover behavior from quasi-equilibrium 
to off-equilibrium behavior and its scaling properties are, interestingly enough, very similar 
to that found in the spin-autocorrelation function of the 3 dimensional spin-glass model . 
The 2nd is the break down of the fluctuation dissipation theorem (FDT) associated with 
the response to the tilt field. Again the correlation function which is conjugate with the 
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response function shows clear waiting time dependence. On the other hand, in contrast with 
the case of spin-glasses, the response function itself does not show waiting time dependence. 
The FDT which relates the two breaks down at t ^ t^, which coincide with the crossover 
behavior in the dynamical overlap function. 

In our analysis we also utilize the transfer matrix method for the following two purposes. 
Firstly, we check the consistency between the static limit of the data of the dynamical 
quantities and the static expectation values. Secondly we study the relation between the 
complex free energy landscape and the relaxational dynamics. Actually we could visually 
monitor the thermal jumps between the excited states. The structure of the web of the 
excited states is very complex which consists of numerous loop-like structures of various 
sizes (see Fig. 3 below). Since this model is rather simpler than spin-glass models, we 
believe that understanding the dynamics of this model would give valuable insight into the 
glassy dynamics of more complicated systems. 

This paper is organized as the following. In section 2 we describe our model and introduce 
the Monte Carlo dynamics. In section 3, the elementary process of the dynamics is studied 
combining with the analysis of the free energy landscape by the transfer matrix method. 
In section 4, the procedure of the temperature quench experiment is described. The results 
are given in section 5. In section 6 we conclude this paper with some phenomenological 
arguments. 

II. THE MODEL 

We study a lattice version of 1-1-1 dimensional directed polymer in random media 
(DPRM). The Hamiltonian is 

N 

= E {9\<^) -x{z + l)\+ fx{x{z), z)} - hx{N) (1) 

z=l 

where x{z) {z — 1- ■ ■ N) represents the configuration of the polymer whose length is N. As 
shown in Fig. 1, the polymer is situated on a square lattice and directed along the 2;-axis 
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so that overhangs are excluded. We impose the so called restricted sohd on solid (RSOS) 
condition, i.e. the steps \x{z) — x{z — 1)| are only allowed to take integer values —1, and 
1. We fix one end x{0) at x = while we allow the other end x{N) to move freely. The 
1st term in the Hamiltonian represents the elastic energy. The 2nd term represents the 
random potential which takes random numbers distributed uniformly over — u < < u 
independently on every lattice site on the square lattice. The last term is the tilt field which 
drives the free end and tilts the polymer. 

We model the relaxational dynamics of the present model by an ordinary heatbath Monte 
Carlo dynamics as in the following way. In one micro Monte Carlo step (micro-MCS), one 
segment, say the segment at z — Zq is tried to move. The new position is chosen among the 
possible choices allowed by the RSOS condition with the appropriate transition probabilities 
so as to satisfy the detailed balance condition. As shown in Fig. 2, this micro-MCS can be 
classified into three cases depending on the positions of the neighboring segments at — 1 
and Zq + I, namely case a), b) and c) which corresponds to \x{zo — 1) —x{zo + 1)| =0 ,1 and 
2 respectively. The entire configuration is refreshed by vectorized sublattice fiippings in one 
Monte Cairo step (MCS) (=N micro-MCS). One sublattice consists of the even numbered 
segments z — 2,4,6... and the other sublattice consists of the odd numbered segments 
2; = 1, 3, 5 . . .. 

III. ELEMENTARY PROCESSES OF THE DYNAMICS 

In this section we consider the elementary processes of the relaxational dynamics in the 
present model. In any quenched random systems, it is generally essential to understand the 
structure of the free energy landscape in order to study the slow relaxational dynamics, which 
is unfortunately very hard to do on spin-glass models. However, thanks to the simplicity 
of the model, we can try such an attempt on the present model in the following way. We 
consider the spatial variation of Pf^iz, x): the probability to find the polymer going through 
a particular lattice point {z, x) on a certain sample of random potential /i. It is defined as 
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P^{z,x) = Z-^ (5(x(^)-x)(5(x(0)-0)exp(-/37^[x,/i]) (2) 

configuration 

where the sum is taken over the all possible configurations (or paths) {a;(l), . . . ,x(A^)} and 
is the normalization factor defined so that I]x=-2 -^^(-^5 ~ ^- straightforward to 
calculate this probability by transfer matrix method In Fig. 3, an example of the spatial 
variation of P^{z^x) of a system of = 100, T = 0.3, g = 0.01 and a = 0.5 on a certain 
sample of random potential /x is displayed by a density plot; the intensity increases as the 
color becomes brighter. We recognize that most of the probability density is confined in the 
white tubes which are understood as the thermally active excited states. In the dynamics, 



they presumably serve as traps |T^, which tend to trap the polymer for long times. On 
the other hand there are black voids between the tubes which presumably serve as free 
energy barriers. It is quite remarkable that the tubes form a very complicated network 
which consists of various sizes of loop-like structures. It is qualitatively consistent with the 
prediction by Hwa and Fisher Q who have predicted a broad distribution of the size of the 
loops. As the temperature increases, one can observe that the diameter of the tubes grows 
and the structure itself changes on large scales as predicted in 0. 

From the above observation, we speculate that the dynamics of the polymer roughly 
consists of the two kinds of elementary processes. One is the fast fiuctuation within each 
tubes (traps). The other is the thermally activated jumps between such tubes, which is 
analogous to the droplet excitations in spin-glasses. Actually these two different processes 
can be visually monitored rather easily as the following. In Fig. 4 we plot the time sequence 
of the segment a;(80) on the same system shown in Fig. 3. The data are averaged over every 
interval of 10^ (MCS) trying to mask small scale fluctuations. Apart from the remanent small 
scale fluctuations, which is understood as the fast fluctuations within each traps (tubes), 
the coarsegrained process is apparently understood as jumps between the four traps which 
correspond to the four major peaks of P{80,x) (white tubes in Fig. 3 at 2; = 80). 

The thermal jumps between the traps must go over the free energy barriers which lie 
between them (black void in Fig. 3). The scaling property of the relaxation times of the 



5 



thermal jumps between the tubes has already been studied in the context of the aforemen- 
tioned transport problem in the vortex glass phase 0. For the convenience of the later 
discussions, we summarize it here. For example, suppose that there is a loop-like structure 
of tubes as schematically shown in Fig. 5. Let its transverse size A. Then the free energy 
barrier -Fb(A) associated with the thermal excitation over the loop, which we hereafter call 
loop- excitation iollowing P|, scales as 

Fb(A) ^ a(T)A", (3) 

where a(T) is some temperature dependent constant. Note that near the free end of the 
polymer, the tubes does not form complete loop structures but rather ?7-shaped structures 
in which one side of the ends of the tubes are left open. But we expect that such thermal 
excitations scales in the same way as the complete loop excitations. The value of the exponent 
a is expected to be 1/2 in 1+1 DPRM Then the relaxation time of the loop-excitation 
scales as 

r(A) ~ exp(FB(A)/T) ^ exp( A^). (4) 

The broad distribution of the size of the loops Q implies broad distribution of relaxation 
times. 



IV. PROCEDURE OF THE TEMPERATURE QUENCH EXPERIMENT 

The numerical experiment is done in the following way (see also the illustration in Fig. 6). 
At first, a random initial configuration is prepared. Then the system is let to evolve for 
(MCS) by the heatbath Monte Carlo dynamics of a certain temperature T. This means that 
the polymer is forced to approach the equilibrium state of temperature T being attached to 
a heat bath of temperature T. After this t^ (MCS) of aging, the polymer's configuration is 
stored and copied to a replica system. Then we let the two systems, say replica A and B, 
continue the relaxational dynamics using the same random numbers (same thermal noise of 
the heat bath) but appling small tilt field h to system B. 
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We measure the following dynamical overlap function as a probe to the relaxational 
dynamics of the unperturbed system (replica A). It is defined as 



q{t^ + t,t^) = /^l^Y.^6{x(^){z,t^ + t),x(^){z,t^))^, (5) 

where 5{x, y) is the Kroneker's delta and x^'^^z^ t) [z = 1 . . . N) represents the configurations 
of the polymer at t (MCS). Here the bracket (. . .) denotes the avarage over different sequences 
of the random numbers and the over bar the average over different realizations of the 
random potentials. It is defined in an analogy with the spin-autocorrelation function of 
the spin-glass models and conveniently measure the 'closeness' between the configurations 
at two different times + t and t„. In the static limit, it is expected to converge to the 
expectation value of the overlap grep between two replica systems which has been introduced 
by Mezard U, 

?rep = Astatic = Jim lim q{t^ + t, t^). (6) 

I— >00 tw— >00 

We measure the dynamical response to the tilt field at t (MCS) after t„ (MCS) of aging 
as the distance between the temporal positions of the free end x{N) of the two replicas A 
and B, which we denote as x['^''{t„ + t) and x[^\t^ + t), 



6xe{t; U) = (xi^\u + t)- xi'^\u + t)). (7) 

Note that (. . .) encloses both x^^^\t) and x^^\t) because we use the same random numbers in 
replica A and B in our procedure. This definition is slightly different from the conventional 



one (^xi^\t^ -|- t)^ — (^xi'^\t^ + ^))- However our present choice is more efficient in practice 
because it provides essentially the same result but with less noise. 

In equilibrium, the response function (|^) is related with the correlation function 



Ceituh) = (xi''\h)xi^\t2)) (8) 



as 



6xe{t; K)=h R^{t + t^, t)dt = f3h [Ce{t^ + t, t„ + t) - C^{t^ + t, t^)] , (9) 



by the fluctuation disspation theorem (FDT) of the 1st kind, 

Reih, t2) = 6x,{U-t2]t2)/6h{t2) = (3dC / (10) 

provided that h is small enough (linear response regime). Note that in the static limit 
^ oo and t ^ oo of (^, we get the ordinary static FDT, 

(xe). = lim lim Sxe{U;t) = I3hx (11) 

t — >00 tw — ^OO 

where the susceptibihty x is given by 

X = {xD - {x,f. (12) 

We performed a series of temperature quench experiments on the system with various 
sizes = 20 . . . 300, and temperatures at T = 0.1 . . . 0.5. The parameter g and a are flxed 
at g = 0.01 and a = 0.5. The average over 10^ . . . 10^ samples of random potential were done 
depending on the system size. In order to ensure that h is small enough so that the linear 
response condition is satisfled, we calculated the disorder averaged static susceptibility x (see 
(0)) and (xe)^ by the transfer matrix method and checked that ( p^Tj ) is well satisfled. 

V. RESULT 

A. Crossover Behavior of Dynamical Overlap Function 

We flrst present the results of the dynamical overlap function deflned in (^ on the 
unperturbed system (replica A). The following analysis is done on the data of = 300 at 
0.1 < T < 0.5. The observation was done in the time window (0 < t < 10"^ (MCS)) for 
i„ = 2, 2^, . . . , 2^° (MCS). The data at T = 0.1 and 0.3 are plotted by double logarithmic 
plot in Fig. 7. The remarkable feature is that the data curves show strong waiting time 

dependence and each of them drops off rapidly at around t It is presumably 

understood as a manifestation of the crossover from the quasi-equilibrium behavior {t <^ t^) 
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to the off-equilibrium behavior (t S> tw), i-e. aging effect. Actually we see later that it 
coincides with the FDT breaking which also occurs at around t ~ t^- 

The quasi-equilibrium regime {t <^ tw) begins with an initial fast decay and crossovers 
to a much slower decay, which is only visible on the data curves of sufficiently large t^. 
Unfortunately we cannot determine the asymptotic functional form {t^ ^ 1) of the latter 
slow decay at this moment. However the scaling analysis which we discuss later suggests 
an algebraic law (see ([T7|) below) with very small exponent x{T) which varies with the 
temperature as shown in Fig. 10. On the other hand, the off-equilibrium part is well fitted 
by an algebraic law 

g(t^,t^ + t) ~r^(^) (t>tw), (13) 

where the exponent A(T) also depends on temperature as shown in Fig. 10. It is interesting to 
note that these characteristics of the crossover behavior are very similar to those originally 
found by Rieger in the spin-autocorrelation function of the 3 dimensional Edwards- 
Anderson spin-glass model in the spin-glass phase. 

We now analize the scaling properties of the dynamical overlap function. We try the 
scaling form 

g(tw, tw + t) ^ C(tw, T)qT{t/U). (14) 

As an example, the scaling plot of the data of = 300, T = 0.17 is shown in Fig. 8. The 
parameter C(t^,T) was chosen so as to scale the off-equilibrium part (t/t^ ^ 1) as good 
as possible. It turns out that the initial fast decay part, which we mentioned before, does 
not fit on the master curve of the t/tw scaling. In practice, we could obtain the scaling plot 
shown in Fig. 8 discarding the data of t < 10 (MCS). At higher temperatures we have to 
discard more data to obtain good master curves. 

As expected, the scaling function behaves as qriy) ~ y^^^'^^ aX y ^ 1 with A(T) which 
we found in (pTSp. On the other hand the constant C(t„, T) turns out to be well fitted by an 
algebraic law of tw, 
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C{U,T)^tl^^\ (15) 

as shown in Fig. 9. The exponent x{T) obtained by this fitting varies with temperature as 
presented in Fig. 10. At T < 0.17, x{T) becomes too small to be determined precisely. 
Using (p!5D, we can rewrite ([T^) as 

+ ~r^(^)0T(tAw) (16) 

where (priv) is related with qriy) in ( P^ by (priy) = y^^^^Qriy)- Note that this scaling form 
is also identical to that originally found in 3 dimensional spin-glass model by Rieger |Q. The 
new scahng form (0) implies that the quasi-equilibrium part (t/t^ ^ 1) behaves as 

qiU, t + tw) ~ t""^^^ (tw > 1 , t/U < 1). (17) 

However as can be seen in the master curve in Fig. 8, the crossover takes place rather 
gradually and the left branch of the master curve (t/tw ^ 1) seem to decay faster than the 
expected algebraic law t~^^'^^ with x{T) obtained from (p^. This discrepancy implies that we 
have not yet attained the asymptotic scaling behavior of the t/t^ part. Unfortunately 
we could not accomplish it with our available computational power. 

Lastly we discuss the finite size effects. As far as the system size N is finite, the size of the 
largest loop-excitation and thus the associated maximum relaxation time Tmax(A^) available 
in the system must be finite (see @)- Thus we expect that system is fully equilibrated if we 
take larger than a certain equilibration time tcq{N) which may be comparable with the 
Tmax(A^)- The finite size effect appears in the dynamical overlap function as the following. In 
Fig. 11 the data on the system of = 20 and T = 0.5 is shown as a typical example. It can 
be seen that the waiting time dependence on the dynamical overlap function saturates at 
some finite as t„ is increased: the curves do not show any further dependence when 

exceeds some characteristic time teq{N). In the example shown in Fig. 11, it is larger 
than 2048 but less than 16384 (MCS). We have checked that teq{N) increases as we increase 
the system size A^. It is also recognized that the dynamical overlap function saturates to 
the static limit ^static (see @) when t ^ tcq- The value of Astatic is nearly equal to q^ep 
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obtained by transfer matrix method which uses two rephcas [^]. Though grep is non-zero in 
finite systems, further calculations by the transfer matrix method on larger systems show 
that grep decrease as the system size increases and it vanish in the thermodynamic limit 
^ oo (see §). 



B. Breakdown of FDT 

We next present the results of the response to the tilt field defined in (|^). The curves in 
Fig. 12 shows the l.h.s and r.h.s. of the FDT relation (^) at different waiting times = 40, 
200 and 1000 (MCS) on the system of = 20 at T = 0.2 with the tih field /i = 0.1 by semi- 
logarithmic plot. The quench experiment was done over 10^ different samples of random 
potential in order to take sufficient average over the disorder. The curves of the l.h.s. of 
(correlation) strongly depends on the waiting time t^. One the other hand, the curves of 
the r.h.s of (||) (response) of corresponding t„ do not seem to have significant dependence 
and almost seem to marge with each other within the numerical accuracy. Comparing with 
the data of larger system sizes, we checked that the data are not spoiled by finite size effects 
within this observation time window (0 < t < 10^ (MCS)). However it is difficult to get 
smooth data on larger systems because they are not self- averaging quantities and sample to 
sample fluctuations are larger on larger systems. 

As is clear from the flgure, the FDT relation is satisfled in the regime t ^ tw but broken 
in the regime t ^ tw It is again a clear evidence of the crossover from quasi-equilibrium 
to off-equilibrium dynamics, i.e. aging effect. Similar FDT breaking was previously found 
in the spin-glass phase of the 3 dimensional spin-glass models by Monte Cairo simulations 
TO| , p!I[] . However there is an important difference: in the spin-glass models both the response 



and correlation functions strongly depends on while in the present model the response 
dose not seem to have dependence. This is understood as the following. As the tilt fleld 
pulls the free end of the polymer, it moves by successive local thermal excitations around it, 
while the rest of the system is indifferent of the existence of the tilt fleld. We expect that 
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the free energy barrier of such thermal excitations scales as (J^). Thus the response will scale 
the same way as the typical transverse size of a loop-excitation which becomes active after 
waiting t, 

5xe(t;tw) ~log(t)^/". (18) 

In order to clarify this scaling form, we performed longer simulation up to 10^ (MCS). In 
Fig. 13 we show the double logarithmic plot of 6xc(t; 0) v.s. log(t) on the system of = 30, 
T = 0.20 and = 60 at T = 0.30 with h = 0.03. It can be seen that the behavior is 
consistent with (|l^) with a = 1/2 and thus supports the above argument. The data on the 
system of different system sizes and at different temperatures show similar behavior except 
that data of smaller systems show saturations within the observation time window. 

Lastly we consider the extremely high temperature case in finite systems. As we men- 
tioned before the diameter of the tubes grow as the temperature increases [Q. Then if the 
system is too small, a single tube swallows up the whole system and the thermal jumps 
between the tubes disappear. In Fig. 14 is shown an example of free energy landscape for 
such case, the density plot of -P^(-2, x) of a system of A^ = 20 at T = 1.0 on a certain random 
potential. Roughly only a single big tube can be recognized. In such an extreme case, we 
found no aging effects. The FDT relation is fulfilled in the whole t range up to the static 
limit as shown in Fig. 15. At around 10^ (MCS), they saturate to the static limit, (xe)^ 
or h/Txe (see (pT|)) calculated by the transfer matrix method. In addition, the dynamical 
overlap function does not show any dependence. These facts suggest that the thermal 
jumps between the tubes, i.e. loop-excitations are the essential processes to bring about the 
observed aging effects. 

VI. CONCLUSION 

The relaxational dynamics of 1+1 dimensional DPRM was investigated by the Monte 
Carlo simulations which mimiced the IRM experiment in spin-glasses. The aging effect 



12 



appears as the systematic waiting time dependence on both the dynamical overlap function 
and the FDT breaking. The dynamical overlap function shows clear crossover from slow 
quasi-equilibrium decay (t <^ t„) to fast off-equilibrium decay {t ^ The latter is well 
fitted by ab algebraic law whose exponent depends on temperature. The t/t^ scaling scheme 
works well and it suggests algebraic decay also in the quasi-equilibrium part. However we 
could not confirm the latter directly with our computational power. On the other hand, 
the response to the tilt field grows as log(t)^ (|18D independently of the waiting time, which 
supports the assumed scaling properties of loop-excitations (see (H) and (H)). Concerning 
the fact that it does not show dependence, this model appears rather trivial compared 
with, for example, spin-glasses in which there are dramatic effects also in the response 
function so that one can actually observe the aging effects in experiments by measuring 
the magnetization None the less, other aspects of the aging effect are very similar 

to that found in the spin-glass phase of 3 dimensional spin-glass models p|- |Tl|] . Hence we 
believe that this model provides a good testing ground for various theoretical ideas of the 
slow dynamics of quenched random systems. 

The important elementary processes to bring about the observed aging effect are the 
loop-excitations, i.e. thermal hoppings between the traps (tubes) in the free energy land- 
scape. They were actually demonstrated by direct monitoring of the time sequences. Thus 
equilibration process is understood as the slow growth of quasi-equilibrium domain driven 
by successive loop-excitations. Note that this growth mechanism of the quasi-equilibrium 
domain is essentially the same as proposed in the spin-glass models in which droplet excita- 
tions drive the growth of the domain 0. From (^), we expect that its transverse size R{t) 
grows as 

R{t) ~ log(t)^/". (19) 

However the scaling properties of the dynamical overlap function are not trivial. Though 
they turned out to be quite similar to those found in the spin-glass model [0, they cannot 
be explained by simply borrowing the assumptions of the above mentioned droplet theory 
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0. In the droplet theory, it is assumed that 

q{t,0) ^ R{t)-\ (20) 

where q{ti, ^2) is now the spin-autocorrelation function and A is some unknown exponent 0. 
(Here A is the notation of and different from A(T) defined in (|TB|).) But we don't have 
any reason to believe that such a relation exists in the present model, and actually our result 
(|13|) cannot be consistent with (^) if we assume (pO|). In the study of spin-glass model 0, 



it is argued that ([T3|) can be consistent with ( pO]) if one assumes that the free energy barrier 
grows logarithmically with its transverse size A rather than algebraically so that R{t) grows 
algebraically with a certain temperature dependent exponent. If it is the also case in the 
present model, the response function (|^) should be fitted by such algebraic law. However 
our data of the response function show systematic curvature in double logarithmic plots 
and do not prefer such algebraic law but the logarithmic law (pISl). Thus we conclude that 
the simple relation (pOf ) does not exist in the present model. The scaling properties of the 
dynamical overlap function require other ways to explain them. 

The important feature we recognize in Fig. 3 is that the tubes are grouped together to 
form some very complicated network. We already know how the relaxation time associated 
with each single loop-excitation scales with its transverse size (see (^). However there is an 
obvious but important rule: a loop-excitation can flip only when some part of the temporal 
configuration of the polymer actually stay in either side of the loop, in other words the 
empty loops can not flip. It should also be noted that the polymer is not allowed to be 
broken into pieces during the dynamical process. Thus not the entire loop-excitations built 
in the system can become active simultaneously but only those which are associated with 
the sequence of the loops which contain the actual temporal configuration of the polymer. 
From the above considerations, it is clear that the information of the connectivity of the 
loops is an important ingredient to describe the dynamics of the polymer. 

Let us now consider the dynamical overlap function (^. It can be roughly interpreted 
as the probability that the segments of the polymer return to the original traps in which 

14 



they initially stay after jumping around other traps in the network. Here we need the 
detailed description of the dynamics mentioned above since the probability is the sum of 
the probabilities associated with all of the possible paths which make such return trips. 
Each of these paths consist of different sequences of loop-excitations which vary in sizes 
and thus have different relaxation times. It is unlikely that there is a single typical mode 
whose probability dominate the total probability. On the other hand, the response to the 
tilt field does not depend directly on such details. The response function scales just as the 
expectation value of the transverse size of the local loop-excitation around the free end which 
can become active after waiting t. It is possible that such expectation value is dominated 
by the contribution from the typical mode whose relaxation time is comparable with t. We 
think this is the reason why we could obtain (|18D using only the scaling property of single 
loop-excitation. 

Although the network looks apparently quite complicated, a most simple characteriza- 
tion of it will be the following. The network structure consists of loops of various sizes which 
are hierarchically nested so that larger ones enclose smaller ones inside. (A similar picture 
was proposed before by Villain concerning the organization of droplets in spin-glass phase, 
see [0). From this simple picture, the dynamics can be regarded as one of the hierarchi- 
cal diffusion processes which have been intensively studied in rather abstract contexts, for 
example the ultradiffusion models [l^. It is interesting to note that in the latter models 
the autocorrelation function (return probability) shows asymptotically algebraic decay with 
some temperature dependent exponent, which is also the case in the dynamical overlap func- 
tion in the present model. Furthermore, it has already been shown by Sibani and Hoffman 



that hierarchical diffusion processes can show aging effects [^. Thus it will be fruitful to 
construct a phenomenological theory of the aging effect of the present model from this point 
of view. A study in this direction will be reported elsewhere . 
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FIGURE CAPTIONS 

Fig. 1: A representation of a configuration of 1+1 dim. RSOS DPRM. 

Fig. 2: The three different ), b) and c) of the local configuration around a segment 

(see text for the definition). The possible new states 1, 2 and 3 of the segment are enclosed 
by boxes. 

Fig. 3: The density plot of P^{z,x) (see text for definition) on a certain sample of random 
potential fi {N = 100, T = 0.3, g = 0.01 and a = 0.5). 

Fig. 4: A time sequence of a;(80) on the same system shown in Fig. 3. The positions of the 
four traps a, b, c and d are indicated. 

Fig. 5: A schematic picture of a loop-excitation. The polymer (bold line) inside the left 
tube jumps over the void to the right tube and vice versa. 

Fig. 6: The procedure of the temperature quench experiment. 

Fig. 7: The waiting time dependence of g(t„, tw + 't)- N = 300, 10^ samples at a) T = 0.1 
and b) T = 0.2. The waiting time increases as = 2,4,8, ...1024 (MCS) from the bottom 
to the top curve. 

Fig. 8: A scaling plot of the dynamical overlap function on the system of = 300 and 
T = 0.17 by the scheme (|14[). The scale of the vertical axis is arbitrary. The doted line 
represents the power law using x obtained by the fit ([15|) . The solid line represents the 
power law t^'^ with A obtained by the direct fit of the t/t^ ^ 1 part. 

Fig. 9: Double logarhithmic plot of C{U, T) v.s. U at T = 0.17, 0.25, 0.30, 0.40 and 0.50. 
The scale of the vertical axis is arbitrary. The dashed lines are the fit by the algebraic law 
(P!3|) with x{T) presented in Fig. 10. 

Fig. 10: Temperature dependence of the exponents : A(T) and x{T). They are obtained 
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from data of (A^ = 300, 10^ samples). 

Fig. 11: Saturation of aging effect due to the finite size effect (A^ = 20, T = 0.5, 500 
samples), grep was calculated by the transfer matrix method. 

Fig. 12: FDT relation (A^ = 20 T = 0.2 10^ samples). The top three curves are the l.h.s. of 
(y) (correlation) of different waiting times t„ = 40,200 and 1000 (MCS). The three curves 
below are the r.h.s of (response) of corresponding t^- 

Fig. 13: Functional form of the response function (A^ = 30, T = 0.20, 10^ samples and 
A^ = 60, T = 0.40, 5 X 10^ samples with h = 0.03). The solid line represents log(t)^ law 

m- 

Fig. 14: The density plot of P{z,x) on a certain sample of small system at a high tempera- 
ture. (A^ = 20, T = 1.00) 

Fig. 15: The FDT relation (^ of a small system at a heigh temperature (A^ = 20, T = 1.00, 
h = 0.1, 500 samples) for t„ = 4, 32, 256 (MCS). The horizontal line is h/Txe where Xe (see 
(0)) is obtained by transfer matrix method. 
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Fig. 1: A representation of a configuration of 1 + 1 dim. RSOS DPRM. 




Fig. 2: The three different cases a), b) and c) of the local configuration around 
a segment (see text for the definition). The possible new states 1,2 and 3 of the 
segment are enclosed by boxes. 
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Fig. 3: The density plot of P^(z,a;) (see text for definition) on a certain sample of 
random potential /i {N = 100, T = 0.3, g = 0.01 and a = 0.5). 
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Fig. 4: A time sequence of a:(80) on the same system shown in Fig. 3. The positions 
of the four traps a, b, c and d are indicated on the right side. 




Fig. 5: A schematic picture of a loop-excitation. The polymer (bold line) inside the 
left tube jumps over the void to the right tube and vice versa. 
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Fig. 6: The procedure of the temperature quench experiment. 
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Fig. 7: The waiting time dependence of q{t^,t^ + t): N 
T = 0.1 and b) T = 0.2. The waiting time increases as 
from the bottom to the top curve. 



300, 10^ samples at a) 
- 2,4,8, ...1024 (MCS) 



1.0 




t " {x=0.034) 

(A,=0.117) 



*10' 10 ' lO" lO' 10^ 10^ lo' 

Fig. 8: A scaling plot of the dynamical overlap function on the system of = 300 
and T = 0.17 by the scheme (14). The scale of the vertical axis is arbitrary. The 
doted line represents the power law using x obtained by the fit (15). The solid 

line represents the power law t^^ with A obtained by the direct fit of the tjty, ^ 1 
part. 
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Fig. 9: Double logarhithmic plot of C{t^,T) v.s. at T = 0.17,0.25,0.30,0.40 
and 0.50. The scale of the vertical axis is arbitrary. The dashed lines are the fit by 
the algebraic law (15) with x{T) presented in Fig. 10. 
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Fig. 10: Temperature dependence of the exponents: A(T) and a;(T). They are 
obtained from data of (iV = 300, 10^ samples). 
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Fig. 11: Saturation of aging effect due to the finite size effect (TV = 20, T = 0.5, 500 
samples), q^^^ was calculated by the transfer matrix method. 
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Fig. 12: FDT relation (iV = 20 T = 0.2 lO'' samples). The top three curves are the 
l.h.s. of (9) (correlation) of different waiting times U = 40,200 and 1000 (MCS). 
The three curves below are the r.h.s of (9) (response) of corresponding ty,. 
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Fig. 13: Functional form of the response function (N = 30, T = 0.20, 10^ samples 
and iV = 60, T = 0.40, 5 x 10^ samples with h = 0.03). The solid line represents 
\og{tf law (18). 
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Fig. 14: The density plot of P{z, x) on a certain sample of small system at a high 
temperature. [N = 20, T = 1.00) 




Fig. 15: The FDT relation (9) of a small system at a heigh temperature {N = 20, 
T = 1.00, h = 0.1, 500 samples) for = 4,32,256 (MCS). The horizontal line is 
h/Txe where Xe (see (12)) is obtained by transfer matrix method. 



